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^ ■ Abstract 

" ' Fermionic path integral Monte Carlo simulations have been applied to study the equilibrium properties of the 

hydrogen and deuterium in the density and temperature range of 1.6 < Va < 14.0 and SOOOiiT < T < 167000ii'. 

C/3 ' We use this technique to determine the phase diagram by identifying the plasma, the molecular, atomic and 

Cy I metallic regime. We explain how one can identify the phases in the path integral formalism and discuss the 

f-«) • state of hydrogen for 5 points in the temperature-density plane. Further we will provide arguments for the 

• ' nature of the transitions between the regimes. 

o 

c/2 , 1 Introduction 

>>■ 

The phase diagram of hydrogen has been studied intensively with different theoretical ap- 
proaches |jll,|Q], simulation techniques [^,0] and experiments ||5|,[^. Prom theory, the prin- 
cipal effects at low densities are well-known. On the other hand, the properties at inter- 
mediate density are not yet well understood, and the phase diagram is not yet accurately 
determined. In particular, the nature of the transition to a metallic state is still an open 
question. 

r^ [ In this article, we would like to show how these questions can be addressed by path 

^D ' integral Monte Carlo (PIMC) simulations. Using this approach, we derived the phase dia- 

gram in Figjl| where we distinguish between molecular, atomic, metallic and plasma regimes. 
We will demonstrate how these different states can be identified from PIMC simulations. 
The imaginary-time path integral formalism |Q] is based on the position-space density ma- 
trix p{TZ,TZ' , f3), which can be used to determine the equilibrium expectation value of any 
"^ . operator O, 

>>: /n\ _T^Op _ fdndn'p{n,n',(3) {n\d\n') 

%;. ^ ' TVp Jdnp{n,n,p) ^ > 

where 7^ represents the coordinates of all particles. The low temperature density matrix 
p{TZ, TZ' ,(i) = {TZ\e~^'^\TZ') can be expressed as product of high temperature density matrices 
p{TZ, TZ, t) with the time step r = j3/M. In position space, this is a convolution, 

p{TZo,'R.m;P) = ■■■ dRidn2 ■■■ dUu-i p(7^o,7^l;r) p(7^l,7^2;'^) ••• p{TZm-i,TZm;t). 

(2) 
This high dimensional integral can be integrated using Monte Carlo methods. Each particle 
is represented by a closed path in imaginary time. Fermi statistics is taken into account by 
considering the fermion density matrix, which can be expressed by considering all permu- 
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Figure 1: The computed phase diagram of deuterium is shown in the temperature-density 
plane. (x,o, A) indicate our PIMC simulations and distinguish between different degrees of 
degeneracy of the electrons (x less than 10% exchanges, o more 10% and A over 80%). The 
four main regimes, molecular, atomic and metallic fluid as well as the plasma are shown. 
The thick solid line specifles the plasma phase transition predicted in [p. The thin solid 
lines specify the approximate location of isobars. 



tations V of identical particles, 



(3) 



V 



where A is the antisymmetrization projection operator. Cancellation of positive and neg- 
ative contributions leads to the ferniion sign problem, which is solved approximately by 
restricting the paths within a nodal surface derived from the free-particle density matrix 



2 Phase diagram of hydrogen and deuterium 

We used PIMC simulation with 32 protons and 32 electrons and a time step r = 1/10^ K 
to generate the phase diagram shown in Fig. ffl. In the low density and low temperature 
regime, we find a molecular fluid. In the proton-proton correlation function shown in Fig. ^, 
one flnds a clear peak at the bond length of 0.75 A. We determine the number of molecules 
as well as other compound particles by a cluster analysis based on the distances. Using 
this approach we can estimate the number of bound states (see [|lO[)- We can also estimate 
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Figure 2: Distribution functions for a selection of 5 simulation of hydrogen at different tem- 
peratures and densities, one in each column: 1) a plasma, 2) a atomic fluid, 3) a molecular 
fluid, 4) a molecular fluid with metallic properties, and 5) a metallic fluid. The rows show 
the following: (1) proton-proton correlation function gpp{r) multiplied by the density, which 
means the area under the peak at the bond length of r = 0.75A indicates the number of 
molecules, (2) [gpeir) — l]r^, where the first peak hints to the existence of bound electrons in 
the ground state, (3) pair correlation function for electrons with parallel spins demonstrat- 
ing the Pauli exclusion principle, (4) pair correlation function for electrons with anti-parallel 
spins, where the peak is caused by a localization of wave function along the molecular bond, 
and (5) distribution of the fraction z^ of electrons involved in a permutation. A peak near 
1^ = represent a small degree of degeneracy of the electrons, while one near v = 1 implies 
a highly degenerate electron gas. 



104 Contrib. Plasma Phys. 39 (1999) 

the fraction of molecules and atoms to determine the regime boundaries. However at high 
density, a clear definition of those species is difficult to give. 

Starting in the molecular regime, one finds that increasing temperature at constant 
density leads to gradual dissociation of molecules followed by a regime, with a majority of 
atoms. The atoms are then gradually ionized at even higher temperatures. Lowering the 
density at constant temperature leads to a decrease in the number of molecules, or atoms 
respectively, due to entropy effects. 

If the density is increased at constant temperature, pressure dissociation diminishes 
the molecular fraction. This transition was described by Magro et. al. [y. Its precise 
nature is still a topic of our current research. Using PIMC simulations, one finds it occurs 
within a small density interval and we predict that it is connected with both the molecular- 
atomic and insulator-metal transition. We determine the fraction of electrons involved in a 
permutation as an indication of electronic delocalization. Permuting electron are required 
to form a Fermi surface, which means that a high number of permutations indicate a high 
degree of degeneracy of the electrons. Permuting electrons form long chains of paths and 
therefore occupy delocalized states. This delocalization destabilizes the hydrogen molecules. 
Before all bonds are broken, one finds a molecular fluid with some permuting electrons, which 
could indicate the existence of a molecular fluid with metallic properties. 

The boundaries of the metallic regime are determined by two effects. With increasing 
temperature, the degree of degeneracy of the electrons is simply reduced. If the temperature 
is lowered, the attraction to the protons becomes more relevant, which localizes the electron 
wave function and decreases the degree of degeneracy also (see Fig. [l|) . 
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